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We investigate the percolation properties of equatorial strips in the two-dimensional 0(3) nonlinear a model. 
We find convincing evidence that such strips do not percolate at low temperatures, provided they are sufficiently 
narrow. Rigorous arguments show that this implies the vanishing of the mass gap at low temperature and the 
absence of asymptotic freedom in the massive continuum limit. We also give an intuitive explanation of the 
transition to a massless phase and, based on it, an estimate of the transition temperature. 



1. Introduction 

This talk, though scheduled in the session on 
Perturbation Theory (PT), does not deal with PT 
as such. But one of its conclusions is that in cer- 
tain nonabelian models the perturbative Callan- 
Symanzik /3-function gives the wrong picture of 
the renormalization group flow. 

In 1991 we developed a rigorous criterion for 
the existence of a massless phase in 2D spin mod- 
els based on percolation properties of certain sub- 
sets of the target spin manifold jjj]. This was 
presented at the 1992 lattice conference together 
with nonrigorous arguments that led to the con- 
clusion that all 2D O(N) most likely had a soft 
low temperature phase, contrary to prevailing ex- 
pectations. 

Advances in computer technology have made 
it feasible to tackle directly the question whether 
the percolation properties leading to the existence 
of a massless phase hold or not. In this talk I re- 
port our results giving direct numerical evidence 
for our old conjecture of the existence of a mass- 
less phase in models for which PT predicts mass 
generation and asymptotic freedom (AF); these 
results were first presented in 

For concretenss and simplicity we are dealing 
with the 2D 0(3) model, a.k.a. classical Heisen- 
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berg model. Its standard version is defined on 
the square lattice, with the standard Hamiltonian 
(action) 

H=-^s(x)-s(y) (1) 

where s(.) is a classical spin (unit 3- vector) liv- 
ing on the vertices x of the lattice and the sum 
is over nearest neighbors (xy); the spins are dis- 
tributed according to the Gibbs measure with 
density exp(—f3H). For the purpose of this in- 
vestigation we modify this model by introducing 
a constraint limiting the size of the deviation al- 
lowed between neighboring spins: 

s(x) ■ s(y) > c (2) 

with — 1 < c < 1 (the so-called cut action with 
cut c), and we also replace the square by a tri- 
angular lattice; the Gibbs measure thus contains 
^-functions enforcing the constraint (||). These 
modifications are made for purely technical rea- 
sons and we have checked that they do not change 
the universality class j|. 

2. The percolation criterion 

We briefly sketch the percolation criterion de- 
veloped in 1991 0: We divide the sphere S2 into 
three pieces: 
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• 'equatorial strip' 5 £ , denned by \s-n\ < e/2 
for some fixed unit vector n. 

• 'upper polar cap' P+, defined by s-n > e/2, 

• 'lower polar cap' V~ , defined by s ■ n < 
-e/2. 

and denote the corresponding subsets of the lat- 
tice for a given configuration by Sf 1 etc.. These 
subsets fall into connected components called 
'clusters'; their mean size we denote by (S e ) etc. 
If this mean size is finite, we say that the clusters 
'form islands'; if there is an infinite cluster, we 
say the subset percolates. There is a third possi- 
bility: that the clusters have divergent mean size, 
but none of them is infinite; this we call 'forma- 
tion of rings'. 

The main result of [jl) is the following 

Theorem: If for a certain c > 1 — e/2 S e does 
not percolate, the model has no mass gap. 

The reason why this is so can be understood by 
using the Fortuin-Kasteleyn (FK) representation 
Q for the imbedded Ising spins a x = sgn(s(x)) 
(which is also the basis of the cluster algorithm) . 

But first one has to notice that in 2D it is not 
possible that two disjoint clusters both percolate, 
and therefore, if c > 1 — e 2 /2, the union of the 
polar caps cannot percolate, because then both 
of them would percolate. But if also S e does not 
percolate, as assumed, a lemma of Russo as- 
sures us that the clusters of each of the three sets 
have divergent mean size. 

But due to the inequality c > 1 — e 2 /2, each 
of the clusters of P+ or P~ has to be contained 
in its entirety inside a FK cluster, and hance also 
the FK clusters have divergent mean size. 

Since the mean size of the FK clusters is equal 
to the susceptibility of the imbedded Ising spins, 
this divergent mean size is incompatible with ex- 
ponential clustering, and thus there is no mass 
gap. 

3. Percolation properties: numerical study 

We investigated numerically the percolation 
properties of the equatorial strip for the special 
case (3 = 0. In this case the parameter c replaces 
the temperature in determining how ordered or 
disordered the system is. 



The results of our investigation j2j are sum- 
marized in the 'percolation phase diagram' Fig.|l|. 
The diagram is semiquantitative, but qualita- 
tively correct, as we will explain. 
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Figure 1. Phase diagram of the 0(3) model on 
the T lattice 



In this figure the solid line is the curve c = 
1 — e 2 /2; above that line the two polar caps can- 
not 'touch' and therefore their union cannot per- 
colate. The dashed line separates a regime in 
which the strip S e percolates (above) and one 
in which it does not (below). For small c, the 
strip forms islands for e below that line. Around 
c = 0.4 a dotted line branches off; below it the 
strip still forms islands, whereas between the dot- 
ted and the dashed lines the clusters of the strip 
have mean infinite size without percolating ('for- 
mation of rings'). The interesting region is the 
one between the solid and the dashed lines: here 
the strip does not percolate but the inequality of 
the theorem holds. So in this regime our theorem 
can be applied and allows us to conclude that for 
c > Co, c Q w 0.7 there is no mass gap. 

Fig|l] also shows that for e < e Q , with e Q w 0.76 
the equatorial strip does not percolate for any c. 
[Alles et al || published a study showing that for 
e = 1.05 and (3 = 2.0 (standard action) the equa- 
torial strip percolates. This is correct, but since 
their choice of parameters is such that they are 
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Figure 2. Ratio (P e )/(S e ) for various e values 
versus c 



both in the massive phase and in the percolation 
region of the strip, it is not very relevant for our 
problem] . 

Let me now explain from which facts our 'phase 
diagram' was derived: In Fig.|] we show that ra- 
tios r = (P e )/(S e ) as a function of c for several 
values of e between .78 and .91 and for lattice 
sizes L from 160 to 1280. 
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Figure 3. Ratio (P .5>/<5 .75> 



It is seen that for small c (depending on e) r 



increases sharply with L (we have data which are 
off the scale of this figure and show that the in- 
crease continues). This expresses the fact that in 
this regime S £ percolates, wheras its complement 
P c forms islands of finite size. At a certain value 
of c the curves for different L intersect and r be- 
comes scale invariant; this is the critical point of 
percolation for the chosen value of e in which both 
sets form rings. The pair (c, e) defines a point on 
the dashed curve in Fig.|l|. 

If we increase c beyond the intersection point, 
the size dependence of r is reversed, indicating 
that we are now in the regime of percolation of 
P. 
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Figure 4. ln(Po.9s) vs InL for c = and c = — 1 



Increasing c still further, the curves come to- 
gether again and remain together, indicating that 
for all c in that regime (depending on e) we are 
in the regime of rings formation of both sets, i.e. 
we have entwred the regime between the dashed 
and dotted curves in Fig. |]. For c — > 1 the ratio 
r converges to the geometric ratio (2 — e)/e of the 
sets P e and S £ . 

To corroborate that for e less than about 0.76 
S e does not percolate, no matter what c is, we 
also measured the ratio r' = (P e t}/(5 , e ), where 
e' is chosen in such a way that the two sets have 
equal area, i.e. e' = 2 — 2e = 0.5. Fig.|| shows that 
for c < 0.4 r' increases sharply with L. Since the 



4 



polar cap cannot percolate, this means that P e ' 
forms rings of arbitrary size, whereas S e forms 
islands of finite size. The behavior changes dras- 
tically around c = 0.4, the L dependence of the 
ratio r' becomes a much milder increase, compat- 
ible with a power law behavior of both numerator 
and denominator, but clearly ruling out percola- 
tion of S e . The only possible interpretation is 
that both sets form rings of arbitrary size. 

The crucial fact for our conclusion is obviously 
that a polar cap can form rings of arbitrary size 
even though it is smaller than a hemisphere. This 
ring formation then prevents percolation of the 
corresponding equatorial strip, and thus allows 
the application of our theorem. 

Our approach has been by necessity in the 
spirit of finite size scaling, studying how various 
quantities change with increasing size. So is it 
conceivable that we arc deceived by finite size be- 
havior that changes its character at some astro- 
nomical lattice size? Obviously in the truly in- 
teresting region near the critical point (which is 
around or slightly below c = 0.7) we cannot work 
on lattices of thermodynamic size. But we did 
an additional test at c = 0, where the correlation 
length is about 53 and we can easily go to ther- 
modynamic lattices: We measured directly the 
mean cluster size (P0.95 as a function of L for L 
up to 1280. The results displayed in Fig.|] show 
a linear dependence of ln(Po.95) on ^ n (L), indi- 
cating a powerlike increase of (P0.95) for lattices 
much larger than the correlation length. This is 
in sharp contrast with the behavior at c = — 1, 
also shown in Fig.^, where one can clearly see 
(P0.95) leveling off. 

If we combine this with another fact, which is 
plausible and which we also checked, namely that 
(P e ) is a monotonically increasing function of c, 
we reach the conclusion that this ring formation 
must persist for all c < 1. 



PT in the models showing perturbative asymp- 
totic freedom in various publications and talks at 
lattice conferences (see for instance 0), but our 
percolation study makes the conclusion unavoid- 
able that PT does not give the correct asymptotic 
expansion for the /3-function. 

A question that has been asked is what is 
'driving the transition' to the massless phase. 
This ia also answered in our recent paper 
it can be understood as the transition froma 
gas of instantons to a gas (or liquid) of super- 
instantons, which are the dominant excitations at 
low temperature and create the disorder required 
by the Mermin- Wagner theorem. In Q we used a 
simple-minded energy-entropy argument to give a 
rough estimate of the transition temperature for 
the standard action model; the result {(3 C rt ~ tt), 
is consistent with knwon numerical results. 

REFERENCES 

1. A.Patrascioiu and E.Seiler, J. Stat. Phys. 
69 (1992) 573; Nucl.Phys.(Proc.Suppl), B30 
(1993) 184. 

2. A.Patrascioiu and E.Seiler, Absence of 
Asymptotic Freedom in Non-Abelian Models, 
hep-th/0002153| . 

3. A.Patrascioiu and E.Seiler, paper in prepara- 
tion. 

4. P. W. Kasteleyn and C. M. Fortuin, J. Phys. 
Soc. Jpn. 26 (Suppl.) (1969) 11; C. M. For- 
tuin and P. W. Kasteleyn, Physica 57 (1972) 
536. 

5. L. Russo, Z. Wahrsch. verw. Gebiete, 43 
(1978) 39. 

6. B. Alles, J. J. Alonso, C. Criado and M. Pepe, 
Phys.Rev.Lett. 83 (1999) 3669. 

7. A.Patrascioiu and E.Seiler, Phys. Rev. Lett. 
74 (1995) 1920; Nucl.Phys.(Proc.SuppL) B 
42 (1995) 826. 



4. Concluding remarks 

Our main conclusion is that the 2D 0(3) has 
a transition to a massless phase at low temper- 
ature, contrary to standard lore, which derives 
from the PT calculation of the Callan-Symanzik 
/3-function. We have questioned the validity of 



